Discrete coherent and squeezed states of many-qudit systems 
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We consider the phase space for a system of n identical qudits (each one of dimension d, with d a primer 
number) as a grid of d" x d" points and use the finite field GF(d") to label the corresponding axes. The 
associated displacement operators permit to define s-parametrized quasidistribution functions in this grid, with 
properties analogous to their continuous counterparts. These displacements allow also for the construction of 
finite coherent states, once a fiducial state is fixed. We take this reference as one eigenstate of the discrete 
Fourier transform and study the factorization properties of the resulting coherent states. We extend these ideas 
to include discrete squeezed states, and show their intriguing relation with entangled states between different 
qudits. 

PACS numbers: 03.65.Ta,03.65.Fd,42.50.Dv 
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I. INTRODUCTION 

The concept of phase-space representation of quantum me- 
chanics, introduced in the pioneering works of Weyl Ol, 
Wigner |2], and Moyal is a very useful and enlighten- 
ing approach that sheds light on the coiTespondence between 
quantum and classical worlds. 

Numerous applications of the phase-space methods to phys- 
ical problems have been extensively discussed in the last 
decades [41, |5l |6l |7t]. However, much of this subject is usu- 
ally illustrated in terms of continuous variables, most often 
position and momentum. 

For discrete systems, or qudits in the modern parlance of 
quantum information, things are less straightforward. Since 
the dynamical symmetry group for a qudit is SU(d), one may 
be tempted to interpret its associated phase space as a gener- 
alized Bloch sphere [8, 9], which is supported by the rigorous 
construction of Kostant lllOtl and Kirilov 1 11 ] in terms of coad- 
joint orbits. Even if this picture is quite popular, especially 
when applied to qubits, one can rightly argue that there is a 
lot of information redundancy there and that the phase space 
should be a grid of points, as one could expect for a truly dis- 
crete system. 

Indeed, apart from some noteworthy exceptions | fT2i [Tsl 
[Till , nowadays there is a wide consensus in picturing the phase 
space for a qudit as a d x d grid. This can be traced back 
to the elegant approach proposed by Schwinger lfl5i[T6[[T7ll . 
who clearly recognized that the expansion of arbitrary op- 
erators in terms of certain operator basis was the crucial 
mathematical concept in setting such a grid. These ideas 
have been rediscovered and developed further by several au- 
thors H HH 11^ J2njlljp,j2l, although the contribu- 
tions of Wootters 11251 uoi 1271 12811 and Galetti and cowork- 
ers ||29l|33,llitl are worth stressing. 

To equip this grid with properties analogous to the geome- 
try of an ordinary plane, it turns out essential [32, ,33^ 34] to 
label the axes in terms of the elements of a finite field GF((i) 
with d elements. It is well known that such a field exist only 
when the dimension is a prime or a power of a prime jssll . 
This, of course, gives a special role to qudits in prime dimen- 
sions, but also is ideally suited to deal with systems of n of 



these qudits. 

Once the natural arena is properly established, the next 
question is how to represent states (and operators) on phase 
space. This is done through quasidistribution functions, which 
allow for the calculation of quantum averages in a way that 
exactly parallels classical mechanics. There are, however, im- 
portant differences with respect to a classical description: they 
come from the noncommuting nature of conjugate quantities 
(like position and momentum), which precludes their simul- 
taneous precise measurement and, therefore, imposes a fun- 
damental limit on the accuracy with which we can determine 
a point in phase space. As a distinctive consequence of this, 
there is no unique rule by which we can associate a classical 
phase-space variable to a quantum operator. Therefore, de- 
pending on the operator ordering, various quasidistributions 
can be defined. For continuous variables, the best known are 
the Glauber-Sudarshan P function lf36l [jTlj. the Husimi Q 
function [38[, and the Wigner W function ll39ll . correspond- 
ing to normal, antinormal, and symmetric order, respectively, 
in the associated characteristic functions. In fact, all of them 
are special cases of the s-parametrized quasidistributions in- 
troduced by Cahill and Glauber [40J . 

The problem of generalizing these quasidistributions 
(mainly the Wigner function) to finite systems has also a 
long history. Much of the previous literature has focused 
on spin variables, trying to represent spin states by continu- 
ous functions of angle variables. This idea was initiated by 
Stratonovich [41], Berezin ll42ll . and Agarwal ll43ll . The re- 
sulting Wigner function, naturally related to the SU(2) dy- 
namical_group, has been further studied by a number of au- 
thors 11441 l45i l46l 14711. has been applied to some problems 
in quantum optics 11481 l49ll and extended to more general 
groups [50[. 

However, these Wigner functions are not defined in a dis- 
crete phase space. A detailed review of possible solutions can 
be found in Ref. [51J. Perhaps, the most popular one is due to 
Wootters and coworkers ll52ll53ll54ll . which imposes a struc- 
ture by assigning a quantum state to each line in phase space. 
Any good assignment of quantum states to lines is called a 
"quantum net" and can be used to define a discrete Wigner 
function. In this paper, we show how to introduce a set of s- 
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parametrized functions, in close correspondence with the con- 
tinuous case. We emphasize that, ahhough some interesting 
work has been done in this direction by using a mod d invari- 
ance ifssl l56ll . our approach works quite well for many-qudit 
systems. 

Another essential ingredient in any phase-space description 
is the notion of coherent states fST]. This is firmly established 
for continuous variables and can easily extended for other dy- 
namical symmetry groups [58]. However, for discrete systems 
we have again a big gap waiting to be filled. The reason for 
this can be traced back to the fact that, as brightly pointed out 
in Ref. |59], in the continuum we have one, and only one, har- 
monic oscillator, while in the discrete there are a lot of candi- 
dates for that role, each one surely with its virtues, but surely 
no undisputed champion. 

The strategy we adopt to deal with this problem is to look 
for eigenstates of the discrete Fourier transform ll60ll . For con- 
tinuous variables, they have a very distinguishable behavior 
that is at the basis of the remarkable properties of coherent 
states. We explore this approach, getting a strikingly simple 
family of discrete coherent states (even for many qudits) ful- 
filling all the requirements. 

To put the icing on the cake, we also extend the notion of 
squeezed states for these systems |l6ll] . For a single qudit, the 
resulting states have nice and expected properties. However, 
when they really appear as highly interesting is for multipar- 
tite systems, since they present an intriguing relation with en- 
tanglement. 

In short, the program developed in this paper can be seen as 
a handy toolbox for the phase-space analysis of many-qudit 
systems, which should be of interest to a large interdisci- 
plinary community working in these topics. 

n. PHASE SPACE FOR CONTINUOUS VARIABLES 

In this Section we briefly recall the relevant structures 
needed to set up a phase-space description of Cartesian quan- 
tum mechanics. This will facilitate comparison with the dis- 
crete case later on. For simplicity, we choose one degree of 
freedom only, so the associated phase space is the plane R^. 

The relevant observables are the Hermitian coordinate and 
momentum operators q and p, with canonical commutation 
relation (with h — I throughout) 

[q,p]^li, (2.1) 

SO that they are the generators of the Heisenberg-Weyl alge- 
bra. Ubiquitous and profound, this algebra has become the 
hallmark of noncommutativity in quantum theory. To avoid 
technical problems with the unbounded operators q and p, it 
is convenient to work with their unitary counterparts ll62ll 

U{q) = exp{-iqp) , V{p) = exp{ipq) , (2.2) 

which generate translations in position and momentum, re- 
spectively. The commutation relations are then expressed in 
the Weyl form 

V{p)Lf{q) = e*«P U{q)V{p) ■ (2.3) 



Their infinitesimal form immediately gives (12. Il l, but ( I2.3l l is 
more useful in many instances. 

In terms of tl and V a displacement operator can be intro- 
duced as 

Diq,p)^e-'^P/'Uip)V{q), (2.4) 

which usually is presented in the entangled form D{q,p) = 
exp [i {pq — qp) ] . However, this cannot be done in more general 
situations. 

The D{q,p) form a complete orthonormal set (in the trace 
sense) in the space of operators acting on TC (the Hilbert space 
of square integrable functions on R). The unitarity imposes 
that D^{q,p) = D{-q,-p), and D{0,0) = i. In addition, 
they obey the simple composition law 

b{q^,p^)b{q2,P2) = e^(f b{q^ + 92, <Z2 + P2) ■ 

(2.5) 

The displacement operators constitute a basic element for 
the notion of coherent states. Indeed, if we choose a fixed 
normalized reference state we can define these coherent 
states as Jsl] 

\q,p) ^ b{q,p)\^o) , (2.6) 

SO they are parametrized by phase-space points. These states 
have a number of remarkable properties, inherited from those 
of D{q,p). In particular, D{q,p) transforms any coherent 
state in another coherent state: 

mi,Pi)\q2,P2) ^e'^P''>--'^'P-y^\qi+q2,q2+P2) ■ (2.7) 

We need to determine the fiducial vector \ipo)- The standard 
choice is to take it as the vacuum |0). This has quite a number 
of relevant properties, but the one we want to stress for what 
follows is that |0) is an eigenstate of the Fourier transform (as 
they are all the Fock states) ll63ll . and so is the Gaussian 

Mq) - m = -471 cxp(-gV2) , (2.8) 

in appropriate units. In addition, this wavefunction represents 
a minimum uncertainty state, namely 

(Ag)2 (Ap)2 = i , (2.9) 

where (Aq)^ and (Ap)^ are the corresponding variances. 

Our next task is to map the density matrix g into a function 
defined on R'^. There exists a whole class of these quasidistri- 
bution functions, related to different orderings of q and p. The 
corresponding mappings are generated by an s-ordered kernel 

W^''\q,p) = Tr[gw<^'\q,p)], (2.10) 

where w'-'^'f is the double Fourier transform of the displace- 
ment operator with a weight fixed by the operator ordering 

w^'\q,p) = -TTT-^ I ^M-i{pq' - qp')]D{q' ,p') 

(27r)^ Jr2 

X {'^Jo\b{q',p')\^a)-'dq'dp' , (2.11) 
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and s £ [—1,1]. The mapping is invertible, so that 

^=7^1 ^^''Hq,P)W'^'\q,p)dqdp. (2.12) 
(27r)^ 7r2 

The Hermitian kernels w^*^ (q, p) are also a complete trace- 
orthonormal set and they transform properly under displace- 
ments 

u-(^) {q, p) - b[q, p) u;(^) (0, 0) (g, p) . (2. 13) 

The symmetric ordering (s = 0) corresponds to the Wigner 
function Ty((7,p) and the associated kernel w*^°)(0, 0) is just 
2'P, where V is the parity operator For the antinormal order- 
ing (s — —1), which corresponds to the Husimi Q function, 
w(o)(0,0) = |0)(0|A. 

The quasidistribution functions ( 12.101 ) fulfill all the ba- 
sic properties required for any good probabilistic descrip- 
tion. First, due to the Hermiticity of w'^^\q,p), they are real 
for Hermitian operators. Second, on integrating W'^'^\q,p) 
over one variable, the probability distribution of the con- 
jugate variable is reproduced. And finally, is 
covariant, which means that for the displaced state g' = 
D{qo,Po) gDHqo,Po), one has 

Wl^\q,p) - W^'\q - 9o,P -Po) , (2.14) 

so that these functions follow displacements rigidly without 
changing their form, reflecting the fact that physics should not 
depend on a certain choice of the origin. 

III. SINGLE QUDIT 

A. Discrete phase space 

We consider a system living in a Hilbert space Hd, whose 
dimension d is assumed from now on to be a prime number. 
We choose a computational basis \i) in Hd (^ = 0, . . . ,d — 
1) which we arbitrarily interpret as the "position" basis, with 
periodic boundary conditions \£ + d) = The conjugate 
"momentum" basis can be introduced by means of the discrete 
Fourier transform i64ll . that is 

K~)=^|^), (3.1) 

where 

and we use the notation 

uj{e) e^p{i2Tr£/d) , (3.3) 

uj = exp{i2TT / d) being a dth root of the unity. Whenever we 
do not specify the ranges in a sum, we understand the index 
running all its natural domain. 



Once we have position and momentum basis, the phase 
space turns out to be a periodic d x d grid of points; i.e., the 
torus 'Zd X where is the field of the integer numbers 
modulo d. 

Mimicking the approach in Sec. |ll] we introduce the oper- 
ators U and V, which generate finite translations in position 
and momentum, respectively. In fact, U generates cyclic shifts 
in the position basis, while V is diagonal 

U''\i) = \i + n), = u;{me) ]£) , (3.4) 

where addition and multiplication must be understood mod- 
ulo d. Conversely, U is diagonal in the momentum basis and 
V acts as a shift, which is reflected also by the fact that 

V^TUT'', (3.5) 

much in the spirit of the standard way of looking at comple- 
mentary variables in the infinite-dimensional Hilbert space: 
the position and momentum eigenstates are Fourier transform 
one of the other. Note that the operators U and V are general- 
izations of the Pauli matrices and cr^, so many authors use 
the notation X and Z for them. 
One can directly check the identity 

F™f7" = Lo{mn) J7"t>" , (3.6) 

which is the finite-dimensional version of the Weyl form of the 
commutation relations and show that they obey a generalized 
Clifford algebra 16?]. 

One may be tempted to define discrete position and momen- 
tum operators. A possible way to achieve this is to write ll66ll 

U = exp(-i27rP/d) , V = exp{i2TTQ/d) , (3.7) 

with 

g-E^I^)(^l' ^ = E^~I^K^1- (3.8) 

I i 

However, for finite quantum systems the Heisenberg-Weyl 
group is discrete, there is no Lie algebra (that is, there are no 
infinitesimal displacements) and the role of position and mo- 
mentum is limited. For this reason our formalism is mainly 
based on the operators U and V. 

Next we introduce the displacement operators 

i)(TO,n) =0(TO,n)J7"T>™, (3.9) 

where (j){m, n) is a phase required to avoid plugging extra fac- 
tors when acting with D. The conditions of unitarity and pe- 
riodicity restrict the possible values of 0. One relevant choice 
(for d > 2) that have been analyzed in the literature is L5 Il1 

0(m, n) = a;(2~^ mn) , (3.10) 

where 2~^ is the multiplicative inverse of 2 in Z^. For qubits, 

(/)(m, n) may be taken as (f){m, n) = ±i™". 
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Without entering in technical details, this choice guaran- 
tees all the good properties, in particular the analogous to 
Eq. (lOl i: 

l)(mi, ni)Z)(m2, 712) — uj[2~^{min2 — TO2"i)] 

xZ)(mi + m2, ni + 712) , (3.11) 

and the following relation 

^Y.Dim,n)^V, (3.12) 

where P is the parity operator Pl^) — \ —£) modulo d. Physi- 
cally, this is the basis for translational covariance and this also 
means that D{m, n) translates the standard basis states cycli- 
cally m places in one direction and n places in the orthogonal 
one, as one would expect from a displacement operator. 



B. Coherent states 

Once a proper displacement operator has been settled, the 
coherent states for a single qudit can be defined as 

\m,n) ^ D{m,n)\ipo) , (3.13) 

where 1-00) is again a reference state. These states are also 
labeled by points of the discrete phase space, as it should be. 

A possible choice ll67ll68[| is to use for |?/'o) the ground state 
of the Harper Hamiltonian 116911 

u + W v + v^ 

H = 2 , (3.14) 

2 2 ' ^ 

which is considered as the discrete counterpart of the har- 
monic oscillator with the proper periodicity conditions. While 
such a replacement is interesting, it is by no means unique. 

We prefer to take a different route, pioneered by Galetti and 
coworkers |60]. We use again as a guide the analogy with the 
continuous case and look for eigenstates |/) of the discrete 
Fourier transform, which play the role of Fock states for our 
problem and are determined by 

{£\nf) = ^'m)- (3.15) 

Obviously, the fact that JF* = 1 implies that it has four eigen- 
values: 1, —1, i, and —i. The solutions of this equation were 
fully studied by Mehta fl^ (see also Ruzzi iH). Taking | V'o) 
as the "ground" state (i.e., ^ = 0) one gets 

|^„) = ^^^c.(H)e-5'=^|^), (3.16) 

and the normalization constant C is given by 

C = ^ e-^^"" = ^3 (o|e"^) , (3.17) 



z?3 being the third Jacobi function Note in passing that 
this fiducial state can be alternatively represented as 

The appearance of the Jacobi function in the present context 
can be directly understood by realizing that this function is a 
periodic eigenstate of the discrete Fourier operator with eigen- 
value + 1 and period d. In addition, it plays the role of the 
Gaussian for periodic variables, which makes this approach 
even more appealing |72]. 

We also observe that \ipo) satisfies a "parity" condition: if 
we write it as \ipo) = J^e K)' '^hen q = c_£. This guaran- 
tees that the average values of I] and y in |i/)o) are the same: 

(^ol^/IV'o) = (^o|l>|^o). 

The coherent states ( 13.131 ) have properties fully analogous 
to the standard ones for continuous variables, as one can check 
with little effort. 

The Harper Hamiltonian commutes with the Fourier oper- 
ator \T , H\ = 0. In fact, the state (13.161 ) is an approximate 
eigenstate of ( 13.141 1 in the high-dimensional limit 

i^|0o)^(^-|^ + |;,)|0o), (3.19) 

which provides another argument for its use as a reference. 

Finally, according to the recent results in Refs. iItsIi and 
lf74ll . the following uncertainty relation holds 

(A[/)2 (Ay)2 > J , (3.20) 

where (AC/)^ = 1 — |(0|i/|'0)P [and an analogous expres- 
sion for (AF)^] denotes the circular dispersion, which is the 
natural generalization of variance for unitary operators. One 
can check that |'0o) saturates this inequality, confirming that it 
is also a minimum uncertainty state. 



C. Quasidistribution functions 

The displacement operators lead us to introduce a Hermi- 
tian s-ordered kernel 

w^'^\m,n) = — Lij{nk — ml) D{m, n) 

X (V'o|^Kn)|V'o)-% (3.21) 

which, as w^**-* (g, p) in Eq. ( 12.1 It , appears as a double Fourier 
transform of D with a weight determined by the operator or- 
dering. However, here the parameter s takes only discrete val- 
ues (s = —1, 0, 1). 

These kernels are normalized and covariant under transfor- 
mations of the generalized Pauli group 

m, n) = w'*^ (fc + m, I + n) . (3.22) 
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FIG. 1: Q function ( I3.28t for the reference state \tpo), which plays 
the role of vacuum for continuous states. 

They can be then conveniently represented as 

(to, n) = D{m, n) w'") (0, 0) (m, n) , (3.23) 

where, according to Eq. ( 13.121 ). w*^*) (0, 0) coincides with the 
parity for s — (d ^ 2), as in the continuous case. 

(s) 

The s-ordered quasidistribution functions W~ are gener- 
ated through the mapping 

W^^"\m,n) = Ti[gw'^'\m,n)] , (3.24) 
which is invertible, so that 

e= i Vw"('*Vm,n)VK('*)(m,?i). (3.25) 
a ^-^ 

These functions fulfill all the basic properties required for the 
probabilistic description we are looking for. Let us apply 
them to the reference state |?/;o) (notice that any other coherent 
state is just a displaced copy of this one). The corresponding 
Wigner function (s = 0) can be obtained after some algebra. 
We omit the details and merely quote the final result: 

M^l^„)(m,n) - ^ ^[(2fc-l-2m)n] 

A: p,g£Z 

X exp[-fc + 2TO + gd- (d- l)/2]^ 
X exp[-(7r/d)(fc+pd-d/2)2], (3.26) 

which, in the limit \, can be approximated by the com- 
pact expression 

W^£\(m,n) ^ ||^(-l)'='c.(mfc-nO 

X exp[-7r(fc2 + z2)/(2d)]. (3.27) 
The Q function {s — — V) for the same state It/iq) reduces to 
<3|^o> = l(^o|i^(m,n)|Vo)|', (3.28) 



which exhibits the additional interesting symmetry 

(5|^o)(m, n) = Q|^g)(-n,TO) . (3.29) 

In Fig.[T]we have plotted this Q function for a 31-dimensional 
system. The aspect of the figure confirms the issues one ex- 
pects from a fairly localized Gaussian state. 

IV. MANY QUDITS 
A. Discrete phase space 

Next, we consider a system of n identical qudits, living in 
the Hilbert space Tidn . Instead of natural numbers, it is conve- 
nient to use elements of the finite field GF(d") to label states: 
in this way we can almost directly translate all the properties 
studied before for a single qudit and we can endow the phase- 
space with many of the geometrical properties of the ordinary 
plane ifssll . In the Appendix we briefly summarize the basic 
notions of finite fields needed to proceed. 

We denote by | A) [from here on, Greek letters will label ele- 
ments in the field GF(c?")] an orthonormal basis in the Hilbert 
space of the system. Operationally, the elements of the basis 
can be labeled by powers of a primitive element using, for 
instance, the polynomial or the normal basis. 

The generators of the Pauli group act now as 

C/,|A) = lA + i/), yjA) =x(mA)|A), (4.1) 

where x(A) is an additive character (defined in the Appendix) 
and the Weyl form of the commutation relations reads as 

V^U, = X{[iy) U.V^ ■ (4.2) 
The finite Fourier transform jTSll 

^" A, A' 

allows us to introduce the conjugate basis |A) = -^lA) and 
also we have 

V^=TU^,P^. (4.4) 

In this way, the concepts delineated in the previous section can 
be immediately generalized. For example, the displacement 
operators are 

b{y.,y)^c^{ti,y)U,V^, (4.5) 

where the phase v) must satisfy the conditions 

v) (j>* jy) = 1, (pip., v) (pi-fi, -v) ^ x(-Mi^) , 

(4.6) 

to guarantee the unitarity and orthogonality of D. We also 
impose 0) = 1 and (/)(0, v) = 1, which means that the 
displacements along the "position" axis /i and the "momen- 
tum" axis V are not associated with any phase. 
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For fields of odd characteristics one possible form of this 
phase is 



(4.7) 



and we have then the same composition law as in Eq. ( 13.111 1. 
namely 

X b{fXi+ ^l2,Vi+V2) ■ (4.8) 
B. Coherent states 



Given our previous discussion, it seems reasonable to ex- 
tend the coherent states (|3.13t in the form 



(4.9) 



where |^o) is a reference state to be determined. In the con- 
tinuous case, the extension of coherent states (12.6b to many 
degrees of freedom is straightforward: they are simply ob- 
tained by taking the direct product of single-mode coherent 
states. To reinterpret (14.9b in the same spirit, one needs first to 
map the abstract Hilbert space Tidn, where the n-qudit system 
lives, into n single-qudit Hilbert spaces Tid®- ■ -^Ti-d- This is 
achieved by expanding any field element in a convenient basis 
{(^j} = li • ■ • ; SO that 



'3 "3 I 



where ij E Z^. Then, we can represent the states as 



\£n) = \il,...Jn), 



(4.10) 



(4.11) 



and the coefficients £j play the role of quantum numbers for 
each qudit. 

The use of the selfdual basis is especially advantageous, 
since only then the basic operators (and the Fourier operator) 
factorize in terms of single-qudit analogues 

= J7"i(g). . .®[/"" , Vf, = F™!®- • -(g)!/"" , (4.12) 

and the displacement operators factorize accordingly 

D{fi, v) = i)(TOi, rii) (g) • • • (g) D{m„, n„) , (4.13) 

where rrij , nj G are the coefficients of the expansion of 
/i and v in the basis, respectively. In consequence, the eigen- 
states of the Fourier transform are direct product of single- 
qudit eigenstates and we can write for the reference state 



l*o> -(g)|^o,) 



(4.14) 



where IV^oj) of the form ( 13. 16b for each qudit (with d > 2). 
For qubits, we have Ii76il 



l*o> 



A (|o)+g|i)), 



(4.15) 



with ^ = V2 - 1. 

Unfortunately, the selfdual basis can be constructed only 
if either d is even or both n and d are odd. This means that 
for such a simple system as two qutrits, this privileged basis 
does not exist. Nevertheless, one can always find an almost 
selfdual basis and one can proceed much along the same lines 
with minor modifications (the interested reader can consult the 
comprehensive review Iplll for a full account of these meth- 
ods). 

It is interesting to stress that for n qubits, the reference state 
(14.15b can be elegantly written in terms of the field elements 
GF(2") as follows 



l*n) = 



1 



(l+^2)„/2 



(4.16) 



QeGF(2") 



where the function h{a) counts the number of nonzero coef- 
ficients aj in the expansion of a in the basis. 

The operator transforming from an arbitrary basis {6'^} into 
the selfdual one {dj } is given by 

r= ^ \mi,...,mn){m[,...,m'J, (4.17) 

/ieGF(2") 



where 



(4.18) 



The operator T is always a permutation and plays the role of 
an entangling (nonlocal) operator. 

Let us examine the simple yet illustrative example of a two- 
qubit coherent state. According to Eq. ( 14.161 ). we have 



1 



1 + 



(4.19) 



where cr is a primitive element. The selfdual basis is {a, a^}, 
and we have the representation 



|0) ^ 100) = 



|a) = |10) = 



(4.20) 



\ct') = 101) = 



In consequence. 



= 111) = 



/I 




Vo 




(4.21) 
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In a non selfdual basis, such as {a, a^}, we have 
|0) = |00), W) = \10), |a=^)-|01), lO^Ill) 
and 



(4.22) 



1 




i+e 



The transition operator ( 14.17b turns out to be 



(4.23) 



r 



/O 1 
10 

10 

Vo 1 



(4.24) 



which is nothing but a matrix representation of the CNOT op- 
erator. 



C. Quasidistribution functions 

The displacement operators ( I4.7l l immediately suggest to 
introduce an s-ordered kernel 



(4.25) 



which, in view of Eq. i4.3i . can also be interpreted as a dou- 
ble Fourier transform of D{^, v). We can next introduce s- 
ordered quasidistribution functions through 



(4.26) 



on relation reads as 

= iE^'"'^(^''^)^f (^'^)- (4.27) 



and the inversion relation reads as 

Q 



Due to the factorization of the character in the selfdual ba- 
sis, the kernels w^^^ (/i, v) factorize in this basis 



{li,v) = \[w^'\m,,nj), (4.28) 



and, consequently, also do the corresponding quasidistribu- 
tions 

i 

For the particular case of the Wigner function, one can 
check that 

^ Wg{^l, u) 5^^ati+l3 = XI ^e(/". ^) ^u-a-^t,~l3 , (4-30) 




FIG. 2: Q function for the reference state l^&o), for a system of three 
qutrits. The order in the axes is as follows: a^^, a^"^, ct^^, a, a^, 

21 23 7 15 4 16 6 8 n 9 12 25 24 5 3 

(T^^, a", (T^^, a^°, a'^°, a^*, a"^^, with a the primitive element. 



that is, the sum over a line of slope a is the same as over a line 
of slope —a^^. The sum over the axes /i and ly are thus equal 

X I^g(M, z.) (5,^0 = El^e(M,f^) '5^.0. (4.31) 

Note also, that the Q function reduces to 

Qein, = (fj,, iy\g\fj,, v) ■ (4.32) 

In Fig. |2] we have plotted this Q function for the reference 
state l^o) in a system of three qutrits. The selfdual basis here 
is {cr, <T^,a^} and the primitive element is a solution of the 
irreducible polynomial + 2a;^ + 1 = 0. 



V. SQUEEZED STATES 

Squeezed states constitute a simple nontrivial enlargement 
of the notion of coherent states. In continuous variables, a 
squeezed state is a minimum uncertainty state that my have 
less fluctuations in one quadrature {q or p) than a coherent 
state. They are generated from the vacuum by using the uni- 
tary squeeze operator 



^(s) = exp[— is {qp + pq)] 



(5.1) 



with a subsequent displacement to an arbitrary point in the 
complex plane 



It is easy to check that 

S'(s)gS't(s) = qe' 



Diq,p)S{5) li^o) 



S{s)pS\s 



p e 



(5.2) 



(5.3) 



so that, the operator S{s) attenuates one quadrature and am- 
plifies the canonical one by the same factor determined by the 
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squeeze factor s, which, for simplicity, we have taken as real. 
As a simple consequence of i5.3i one can verify the transfor- 
mations for U (q) and V{p): 

S{b) U{q) S\b) = U^{q) , S{5) V{p) S^b) = V-'ip) ■ 

(5.4) 

For a single qudit, squeezed states have been recently 
considered in detail in Ref. ll6ll] . using an extended Cahill- 
Glauber formalism. Here, we prefer to follow an alternative 
approach and define a squeeze operator as 



(5.5) 



At first sight, this can appear as a rather abstract choice. How- 
ever, notice that 

Sl ?7" Sl ^U"' , V"' = f " , (5.6) 

which is a direct translation of the action ( 15.4b to this discrete 
case. This also means that in the squeezed "vacuum" 



IV'o;s) = Ss\ipo) , 



(5.7) 



the average values of some powers of the displacement oper- 
ators are the same 

{,Po;s\U\^o;s)^{^Po;s\V'"\%;s). (5.8) 

Perhaps, the clearest way to visualize this squeezing is to 
use a quasidistribution, such as, e.g., the Wigner function. If 
Qs = Ss g sl denotes the density operator of a squeezed state, 
we have 



Wg, (to, n) = Wg{sm, s ^ri) 



(5.9) 



whose geometrical interpretation is obvious and is the phase- 
space counterpart of the property (15. 3t . For reasons that will 
become evident soon, we refer to this as "geometrical squeez- 
ing". We also note the following symmetry property of the 
Wigner function 



/ , " Qs 



Wg^ (m, n) 6n,o = ^ Wg^_, (m, n) 6^^ ■ (5.10) 



For many qudits, our developed intuition suggests a direct 
translation of (15.5b in terms of the field elements in GF((i"), 
namely 

5, =5]|A)M|, ^eGF(d"), (5.11) 



in terms of which we can write relations similar to Eqs. (|5.6b - 
(15. 9b . In fact, one can define a squeezed "vacumm" as in 
Eq. (EtJ, i.e., |*o; <r) = "S^l^o)- In Fig. [3] we plot the Wigner 
function for this squeezed state in a system of three qutrits 
with = CT^. 

Nevertheless, now the squeezing acquires a new physical 
perspective: the squeeze operator (15.11b cannot be, in gen- 
eral, factorized into a product of single qudit squeezing op- 
erators. This means that by applying to a factorized state 




FIG. 3: Wigner function for a squeezed "vacuum" state |4'o, ?), for 
a system of three qutrits, with the same order in the axes as in Fig. 2. 



we generate correlations between qudits; i.e., we create entan- 
gled states. The most striking example is of course the n qubit 
case, since there is no single qubit squeezing. 

To understand these correlations consider a general factor- 
ized state 

|*) = ;^Ca|A)= Y1 ...qJ£i,...,£„), (5.12) 

A ce-^,...,ce„ 

and apply (15.11b . The resulting state turns out to be 

5,1^-)= J2 ^^"nei+...+rn„e„Ki,...,C), (5.13) 
ei,--- 

where 

d-l 

j,k=0 



(5.14) 



hjk - tT{0,e,0k) , hk = tr(^-i0fc) 



tr (written in lower case) is the trace operation in the field (see 
the Appendix) and {9j} is the basis. 

As a example let us consider a three-qubit system. Now 
the selfdual basis is {6i = <j^, 62 ~ a^, 6*3 = c^}, where a 
is a primitive element, solution of the irreducible polynomial 

+ X + 1 — 0. The result of applying S^k to the state (15. 12b 
can be expressed in terms of 

•^■^l*) = X! C't^^AlA) = ^ Cp+gCp+rCq\p,q,r) , 
AGGF(23) P,q,r£l,2 

(5.15) 

Sa^l'^) = ^ C„4a|A) = ^ Cp+g+rCp+rCr\p,q,r) . 

In fact, the transformations {S'^.s, ^o-e} generate the same 
entanglement (except for permutations) as S'^.s, while 
{Se2, Sgi} generate the same entanglement (again except for 
permutations) as 5"^. 
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VI. CONCLUDING REMARKS 

In summary, we have provided a handy toolbox for deal- 
ing with many-qudit systems in phase space. The mathemati- 
cal basis of our approach is the use algebraic field extensions 
that produce results in composite dimensions in a manner very 
close to the continuous case. 

Another major advantage of our theory relies on the use of 
the finite Fourier transform and its eigenstates for the defi- 
nition of coherent states. We believe that this makes a clear 
connection with the standard coherent states for continuous 
variables and constitutes an elegant solution to this problem. 
The factorization properties of the resulting coherent states in 
different bases is also an interesting question. 

We have also established a set of important results that 
have allowed us to obtain discrete analogs of squeezed states. 
While for a single qudit, these squeezed states have the prop- 
erties one would expect from our continuous-variable experi- 
ence, for many qudits an amazing relation with entanglement 
appears. 

We think that the techniques presented here are more than a 
mere academic curiosity, for they are immediately applicable 
to a variety of experiments involving qudit systems. 



APPENDIX A: FINITE FIELDS 

In this appendix we briefly recall the minimum background 
needed in this paper The reader interested in more mathe- 
matical details is referred, e.g., to the excellent monograph by 
Lidl and NiedeiTeiter llssll . 

A commutative ring is a nonempty set R furnished with two 
binary operations, called addition and multiplication, such 
that it is an Abelian group with respect the addition, and the 
multiplication is associative. Perhaps, the motivating example 
is the ring of integers Z, with the standard sum and multipli- 
cation. On the other hand, the simplest example of a finite 
ring is the set Z„ of integers modulo n, which has exactly n 
elements. 

A field is a commutative ring with division, that is, such 
that does not equal 1 and all elements of F except have 
a multiplicative inverse (note that and 1 here stand for the 
identity elements for the addition and multiplication, respec- 
tively, which may differ from the familiar real numbers and 
I). Elements of a field form Abelian groups with respect to ad- 
dition and multiplication (in this latter case, the zero element 
is excluded). 

The characteristic of a finite field is the smallest integer d 
such that 



dl = 1 + 1 



1 = 



(Al) 



d times 



and it is always a prime number Any finite field contains a 
prime subfield Zd and has d" elements, where n is a natural 
number Moreover, the finite field containing d" elements is 
unique and is called the Galois field GF((i"). 

Let us denote as [x] the ring of polynomials with coeffi- 
cients in Zrf. Let P{x) be an irreducible polynomial of degree 



n (i.e., one that cannot be factorized over Z^). Then, the quo- 
tient space Zii[X]/ P{x) provides an adequate representation 
of GF((i"). Its elements can be written as polynomials that are 
defined modulo the irreducible polynomial P{x). The multi- 
plicative group of GF((i") is cyclic and its generator is called 
a primitive element of the field. 

As a simple example of a nonprime field, we consider the 
polynomial x^ + x + 1 = 0, which is irreducible in Z2. If 
(T is a root of this polynomial, the elements {0,1,(t, <t^ = 
CT + 1 = (T^^} form the finite field GF (2^) and is a primitive 
element. 

A basic map is the trace 



tr(A) = A + + . . . + A"*" . 
It is always in the prime field and satisfies 
tr(A + A') = tr(A) + tr(A') . 
In terms of it we define the additive characters as 



X(A) = exp 



P 



■tr(A) 



(A2) 



(A3) 



(A4) 



which posses two important properties 

x(A + A')=x(A)x(A'), 



J2 X(AA') = d"<5, 

A'GGF(d") 



'0,A ■ 



(A5) 

Any finite field GF((i") can be also considered as an n- 
dimensional linear vector space. Given a basis {0j}, (j = 
1 , . . . , n) in this vector space, any field element can be repre- 
sented as 



(A6) 



with £j E 'Ld- In this way, we map each element of GF((i") 
onto an ordered set of natural numbers A <^ (£1 , . . . , ^„). 
Two bases {61^ . . . , 9^} and {d[, . . . , 6*,^} are dual when 



(A7) 



A basis that is dual to itself is called selfdual. 

There are several natural bases in GF((i"). One is the poly- 
nomial basis, defined as 



{l,a,a2,...,(T"-i}, 



(A8) 



where cr is a primitive element. An alternative is the normal 
basis, constituted of 



(A9) 



The choice of the appropriate basis depends on the specific 
problem at hand. For example, in GF(2^) the elements 
{cr, cr^} are both roots of the irreducible polynomial. The 
polynomial basis is {1, cr} and its dual is {cr^, 1}, while the 
normal basis {cr, cr^} is selfdual. 

The selfdual basis exists if and only if either d is even or 
both n and d are odd. However for every prime power d^. 
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there exists an almost selfdual basis of GF((i"), which satis- + x + 2 = 0, form an almost selfdual basis 

fies the properties: tr{6i9j) = when i j and tr{9'^) = 1, 

with one possible exception. For instance, in the case of two tr(a^a^) = 1 tr(cr^cr^) = 2 trfcr^cr^) = (AlO) 

qutrits GF(3^), a selfdual basis does not exist and two ele- ' 
ments {a^,a'^}, a being a root of the irreducible polynomial 
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